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Summary. Following on from the recent results of Moradi et al. (2009), we em- 
ploy two established numerical forward models (a 3D ideal MHD solver and MHD 
ray theory), in conjunction with time-distance helioseismology to probe the lateral 
extent of wave-speed perturbations produced in regions of strong, near-surface mag- 
netic fields. We aim to continue our comparative studies of the forward models by 
complimenting our previous surface-focused travel-time measurements through the 
application of a common midpoint deep-focusing scheme that avoids the use of os- 
cillation signals within the sunspot region. The idea here is to also test MHD ray 
theory for possible application to future inverse methods. 



1 Introduction 

In Moradi et al. (2009) we utilized two recently developed numerical MHD 
forward models, in conjunction with surface- focused (i.e., centre-to-annulus) 
time-distance measurements, to produce numerical models of travel-time inho- 
mogeneities in a simulated sunspot atmosphere. The resulting artificial travel- 
time perturbation profiles clearly demonstrated the overwhelming influence 
that MHD physics, as well as phase-speed and frequency filtering, have on 
local helioseismic measurements in the vicinity of sunspots. 

However, there are numerous caveats associated with surface-focused time- 
distance measurements that use oscillation signals within the sunspot region, 
as the use of such oscillation signals is now known to be the primary source 
of most surface effects in sunspot seismology. These surface effects can be 
essentially categorized into two groups. The first revolves around the degree 
to which observations made within the sunspot region are contaminated by 
magnetic effects (e.g., Braun 1997; Lindsey & Braun 2005; Schunker et al. 
2005; Braun & Birch 2006; Couvidat & Rajaguru 2007; Moradi et al. 2009), 
while the second concerns the degree to which atmospheric temperature strat- 
ification in and around regions may affect the absorption line used to make 
measurements of the Doppler velocity (e.g., Rajaguru et al. 2006, 2007). 
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There have been attempts in the past to circumvent such problems by 
adopting a time-distance measurement geometry known as "deep-focusing" 
which avoids the use of data from the central area of the sunspot by only 
cross-correlating the oscillation signal of waves that have a first-skip distance 
larger than the diameter of the sunspot (e.g., Duvall 1995; Braun 1997; Zhao 
& Kosovichev 2006; Rajaguru 2008). In this analysis, we follow up on the 
comparative study presented in Moradi et al. (2009) by using our two estab- 
lished forward models, in conjunction with a deep-focusing scheme known as 
the "common midpoint" (CMP) method to probe the sub-surface dynamics 
of our artificial sunspot. 



2 The Flux Tube and Forward Models 

The background stratification of our model atmosphere is given by an adiabat- 
ically stable, truncated polytrope (Bogdan et al. 1996), smoothly connected 
to an isothermal atmosphere. The truncated polytrope is described by: index 
m = 2.15, reference pressure pq = 1.21 x 10^ g cm~^ s~^ and reference den- 
sity po = 2.78 X 10"''' g cm~^. The flux tube (peak field strength of 3 kG) is 
modeled by an axisymmetric magnetic field geometry based on the Schliiter 
& Temesvary (1958) self-similar solution. The derived MHS sunspot model 
achieves a consistent sound-speed decrease (see Figure 1), with a peak re- 
duction of ^ 45% at the surface {z = 0) and less than 1% at z = — 2 Mm, 
while the one-layered wave-speed enhancement is also confined to the near- 
surface layers, approaching ~ 200% at the surface and around less than 0.5% 
at z = -2 Mm. 
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Fig. 1. Some properties of the model sunspot atmosphere. Left panel: the near- 
surface thermal/sound-speed perturbation profile, shown as a function of sound- 
speed squared. Right panel: shows a Doppler power map, normalized to the quiet 
Sun. 



The two forward models presented in Moradi et al. (2009) are again used 
for our analysis. The first forward model integrates the linearized ideal MHD 
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wave equations according to the recipe of Hanasoge (2008), where waves are 
excited via a prc-computcd deterministic source function that acts on the ver- 
tical momentum equation. To simulate the suppression of granulation related 
wave sources in a sunspot (e.g., Hanasoge et al. 2008), the source activity is 
muted in a circular region of 10 Mm radius. The simulations produce artificial 
line-of-sight (Doppler) velocity data cubes, extracted at a height of 200 km 
above the photosphere, in effect, mimicking Michelson Doppler Imager (MDI) 
Dopplcrgrams. The data cubes have dimensions of 200 x 200 Mm^ x 512 
minutes, with a cadence of 1 minute and a spatial resolution of 0.718 Mm. 
Figure 1 depicts a normalized power map derived from the simulated Doppler 
velocity measurements. 

The second forward model employs the MHD ray tracer of Moradi & Cally 
(2008), where (2D) ray propagation is modeled through solving the governing 
equations of the ray paths derived using the zeroth order cikonal approxima- 
tion and the magneto-acoustic dispersion relation. It should be noted that 
neither forward model accounts for the presence of sub-surface flows. 

3 Common Midpoint Deep-Focusing 

Often utilized in geophysics applications such as multichannel seismic ac- 
quisition (Shearer 1999), the CMP method measures the travel time at the 
point on the surface halfway between the source and receiver (see Figure 2). 
Cross-correlating numerous source-receiver pairs in this manner results in the 
method being mostly sensitive to a small region in the deep interior surround- 
ing the lower turning point of the ray. A re-working of this method has been 
applied to helioseismic observations by Duvall (2003), and has the obvious 
advantage of allowing one to study the wave-speed structure directly beneath 
sunspots without using the oscillation signals inside the perturbed region. 

Our method for measuring time-distance deep-focus travel times is some- 
what similar to the approach undertaken by Braun (1997) and Duvall (2003). 
First, the annulus-to-annulus cross- covariances (e.g., between oscillation sig- 
nals located between two points on the solar surface, a source at ri and a 
receiver at r2, as illustrated in Figure 2) are derived by dividing each annulus 
(A = |r2 — ri|), into two semi-annuli (each being one data pixel wide) and 
cross-correlating the average signals in these two semi-annuli. Then, to further 
increase the signal-to- noise ratio (SNR), we average the cross-covariances over 
3 distances, respectively slightly smaller than, and larger than, A. In the end, 
the 5 (mean) distances chosen (A = 42.95, 49.15, 55.35, 61.65 and 68 Mm 
respectively) are large enough to ensure that we only sample waves with a 
first-skip distance greater than the diameter of the sunspot at the surface 
(~ 40 Mm). 

Due to the oscillation signal at any location being a superposition of a large 

number of waves of different travel distances, the cross-covariances are very 
noisy and need to be phase-speed filtered first in the Fourier domain, using a 
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Fig. 2. An illustration of the CMP dccp-focus geometry indicating the range of rays 
used for this study. The CMP method measures the travel time at the point on the 
surface located at the half-way point between a source (ri) and receiver (r2). For 
the above rays, the CMP is located on the central axis of the spot (r = Mm). 

Gaussian filter for each travel distance. The application of appropriate phase- 
speed filters isolates waves that travel desired skip distances, meaning that 
even though we average over semi-annuli, the primary contribution to the 
cross covariances is from these waves. In addition to the phase-speed filters, 
we also apply an /-mode filter that removes the /-mode ridge completely 
(as it is of no interest to us in this analysis), and we also apply Gaussian 
frequency filters centred at w = 3.5, 4.0 and 5.0 mHz with 6lu = 0.5 mHz 
band- widths, to study frequency dependencies of travel times (e.g., Braun & 
Birch 2008; Moradi ct al. 2009). To extract the required travel times, the cross- 
covariances arc fitted by two Gabor wavelets (Kosovichev & Duvall 1997): one 
for the positive times, one for the negative times. 

Even after significant filtering and averaging, the extracted CMP travel 
times are still inundated with noise. This is certainly an ever-present com- 
plication in local helioscismology as there is a common expectation (with all 
local hclioscismic methods and inversions) of worsening noise and resolution 
with depth. Realization noise associated with stochastic excitation of acous- 
tic waves can significantly impair our ability to analyse the true nature of 
travel-time shifts on the surface (and by extension, also affect our interpre- 
tation of sub-surface structure). But, as we have full control over the wave 
excitation mechanism and source function, we have the luxury of being able 
to apply realization noise subtraction to improve the SNR and obtain statis- 
tically significant travel-time shifts from the deep-focus measurements. This 
is accomplished in the same manner as Hanasoge et al. (2007), i.e., by per- 
forming two separate simulations, one with the perturbation (i.e., the sunspot 
simulation), and another without (i.e., the quiet simulation). We then sub- 
tract the travel times of the quiet data from its perturbed counterpart (see 
Figure 3), allowing us to achieve an excellent SNR. 
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Finally, in order to compare theory with simulations, we also estimate 
deep-focusing time shifts using the MHD ray tracer of Moradi & Cally (2008). 
The single-skip magneto-acoustic rays are propagated from the inner (lower) 
turning point of their trajectories at a prescribed frequency (see e.g.. Figure 
2). These rays do not undergo any additional filtering as the required range of 
horizontal skip distances are simply obtained by altering the depth at which 
the rays are initiated from. The resulting mean (phase) travel-time shifts 
(STmean) derived from both forward models are presented in Figures 3 and 4. 
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Fig. 3. Examples of (phase-speed filtered) CMP mean travel-time perturbation 
(Sr^ean) maps for A = 42.95 (top), A = 49.15 (middle) and A = 61.65 Mm 
(bottom). Left panels: before realization noise subtraction. Right panels: after sub- 
traction. A frequency filter centred at 5.0 mHz has been applied to the data. 
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4 Results and Discussion 

A number of travel time maps derived from the time-distance analysis, both 
before and after noise subtraction, are presented in Figure 3. The impact 
of realization noise subtraction is self-evident in these figures as it is only 
after removing the background noise that we are able to detect statistically 
significant travel-time shifts. The umbral averages of these time shifts are 
shown in Figure 4. The STmean range from a couple of seconds at 3.5 and 
4.0 mHz, to around five seconds at 5.0 mHz. However, even though the size of 
the measured time shifts are significant, there is no clear frequency dependence 
associated with the them. As we arc only using waves outside of the perturbed 
region, surface effects can be effectively ruled out as the cause of the time 
shifts. 

It is worth noting that linear inversions of surface-focused travel time maps 
of actual observations have suggested a two-layered wave-speed structure be- 
low sunspots - a wave-speed decrease of ~ 10 — 15% down to a depth of 
~ 3 — 4 Mm, followed by a wave-speed enhancement, reportedly detected 
down to depths of ~ 17 — 25 Mm below the surface (Kosovichev et al. 2000; 
Couvidat et al. 2006). However, with our forward model clearly prescribing 
both a shallow wave- and sound-speed perturbation profile (Figure 1), it is 
hard to fathom that the time-distance dTmean we are observing can be asso- 
ciated with some kind of anomalous deep sub-surface perturbation. In fact, 
both the sound-speed decrease and wave-speed enhancement at such depths 
registers at less than one-tenth of one percent and the value of the plasma (3 
is in the range of ~ 7 — 18 x 10^ - in all likelihood not significant enough to 
produce a 3-5 second travel-time perturbation. In order to try and identify 
the root cause of these apparent travel-time shifts, it is useful to compare the 
time-distance CMP measurements with those derived from MHD ray theory 
in Figure 4. 

The ray theory CMP 6Tmean clearly appear to be significantly smaller 
at all frequencies, with all observed time shifts registering at less than half 
a second. Certainly, these time shifts are more in line with our expectations 
given the absence of any significant deep sound/wave-speed perturbation. But, 
we must bear in mind the differences between the two forward models before 
drawing our conclusions. With regards to hclioscismic travel times, Bogdan 
(1997) has emphasized that they are not only sensitive to the local velocity 
field along the ray path, but also to conditions in the surrounding medium - 
a clear consequence of wave effects. As such, wave-like behaviour needs to be 
considered when interpreting travel times - something which ray theory does 
not clearly account for, resulting in possible underestimation of deep-focus 
travel times. 

On the other hand, we must also consider the effects of phase-speed fil- 
tering (which is absent in the ray theory calculations) on the time-distance 

measurements. If we look closely at the time-distance STmean maps in Fig- 
ure 3, we notice that they are somewhat smeared in appearance, with the 
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Fig. 4. Observed CMP travel-time shifts as a function of wave/ray travel distance 
(A). Left panel: umbral averages of the CMP time shifts derived from time-distance 
analysis of the simulated data. Right panel: ray theory CMP travel-time shifts de- 
rived from rays propagated at various depths and with a CMP at r = Mm. Light 
sold lines are indicative of frequency filtering centred at 3.5 mHz, dashed lines indi- 
cate 4.0 mHz and bold solid lines indicate 5.0 mHz. 

central sunspot region becoming increasingly sprawled-out across the map as 
we increase A. This behaviour is most likely a consequence of both the phase- 
speed filtering (i.e, the size of the central frequency filter, the filter width, etc., 
see Couvidat & Birch 2006), and the averaging scheme applied to the cross- 
correlations - both of which are a necessity in order to improve the SNR in 
time-distance calculations. These effects, combined with the delocalized na- 
ture of the CMP measurements, may also introduce spurious travel-time shifts. 
However, further testing and control simulations are required to confirm this. 

5 Conclusion 

At the present time, it is sufficient to say that we do not have a definitive 
diagnosis with regards to the above-discussed differences in the size of the 
deep-focus time shifts produced by the two forward models. It may well be 
that we are applying ray theory to regimes where it may be seriously limited. 
On the other hand, the very same could be said about local time-distance 
analysis! Whatever the case may be, these preliminary results have certainly 
provided us with the motivation to conduct further time-distance studies using 
the CMP method. 

The direct (and indirect) effects of phase-speed filtering on deep-focus 
measurements, derived from both simulations and real data, also warrants a 
more detailed examination, as any artifact produced by the filtering process is 
likely to be even more pronounced for phase-speed filtered MDI data, where 
we do not yet have the luxury of realization noise subtraction. These issues are 
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something that we hope to address in the very near future with some ongoing 
comparative studies. 
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